Quasiparticles at the verge of localization near the Mott 
metal-insulator transition in a two-dimensional material 
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The dynamics of charge carriers close to the Mott transition is explored theoretically and exper- 
imentally in the quasi two-dimensional organic charge-transfer salt, k-(BEDT-TTF)2Cu[N(CN)2]- 
Br^Cli-x, with varying Br content. The frequency dependence of the conductivity deviates signif- 
icantly from simple Drude model behavior: there is a strong redistribution of spectral weight as 
the Mott transition is approached and with temperature. The effective mass of the quasiparticles 
increases considerably when coming close to the insulating phase. A dynamical mean-field-theory 
treatment of the relevant Hubbard model gives a good quantitative description of the experimental 
data. 



Understanding strongly correlated electron materials 
such as the high-T c cuprate superconductors, heavy 
fcrmions, transition metal oxides, and organic super- 
conductors is a significant theoretical challenge. Opti- 
cal conductivity measurements are a powerful probe of 
the competition between the itineracy and localization 
of charge carriers 0, 0, H, 0] ■ Unlike in elemental met- 
als, a large redistribution of spectral weight (SW) oc- 
curs with variation in the Coulomb interaction and/or 
temperature [1, 0, H, H, @] ■ The simple Drude model [l| 
fails to describe spectra in these systems, despite its suc- 
cess in conventional metals. In this Letter, we present a 
combined theoretical and experimental analysis of the 
redistribution of SW that occurs as the Mott metal- 
insulator transition is approached in a two-dimensional 
material. 

The K-family of layered organic superconductors based 
on the bis-(ethylenedithio)tetrathiafulvalene (BEDT- 
TTF) molecules [8J have conduction bands which are ef- 
fectively half filled (due to dimerization of neighboring 
BEDT-TTF molecules). As the on-site Coulomb repul- 
sion U is comparable to the bandwidth W, these ma- 
terials are ideal systems to explore the bandwidth con- 
trolled Mott metal-insulator transition (MIT) in two di- 
mensions. Measurements of the DC conductivity have 
been made as the materials tuned through the MIT via 
chemical substitution, pressure, magnetic field, and/or 
temperature @, Much less is known about the 

charge-carriers dynamics in these systems close to the 
MIT @, EH, E EDI- The optical conductivity calcu- 
lated from a dynamical mean- field theory (DMFT) [lil] 
treatment of the relevant Hubbard model is compared 
with optical spectra of single crystals of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br a: Cli_ x . Partial substitution of Br 
by CI in the anion layer increases the effective Coulomb 
repulsion U/W, and drives the system across the metal- 
insulator phase boundary, which is at x ~ 0.7 [HI]. 



The simplest strongly correlated electron model for the 
conducting layers of the k-(BEDT-TTF) 2 X family is a 
Hubbard model on a frustrated square lattice at half fill- 
ing [13: 

i 2(7 

where a is the spin index, t\ describes hopping along one 
diagonal of the square plaquettes, i 2 describes hopping 
around the plaquettes, and cJ CT and c J(T are hole creation 
and destruction operators on the antibonding orbitals of 
BEDT-TTF dimers. The non-interacting model (U = 0) 
has a dispersion relation: 

e(k) = —2t2[cos(k x d)+cos(k y d)]—2tiCOs((k x +ky)d) (2) 

with d the nearest-neighbor dimer distance. We take 
the nearest-neighbor hopping amplitudes to be i 2 = 
—0.03 eV and t\ — 0.8i 2 which leads to a non-interacting 
bandwidth of W ~ 0.3 eV, comparable to values from 
density functional theory (DFT) calculations [lij ]. 

Previous studies of model |T]) based on DMFT cal- 
culations [l|| and its cluster extensions [2(| considered 
the Mott-Hubbard metal-insulator transition driven by 
U/W. Anomalous transport properties observed in k- 
(BEDT-TTF) 2 X [§], such as the non-monotonic tem- 
perature dependence of the DC resistivity, thermopower, 
and Hall coefficient agrees with DMFT, which predicts 
the gradual destruction of quasiparticles with increasing 
temperature, T. This leads to a crossover from Fermi- 
liquid behavior at low T to 'bad' metallic behavior for 
T > T*, with T* « W 0. Signatures of the grad- 
ual destruction of quasiparticles as the temperature is 
increased above T* w 100 K can be seen in Fig. [TJ 
which shows the frequency dependence of the real part 
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FIG. 1: (Color online). Frequency dependence of the opti- 
cal conductivity at several temperatures. Calculations are 
based on a DMFT treatment of the Hubbard model on the 
anisotropic triangular lattice. For U = 0.3 eV, a gradual 
suppression of the Drude peak occurs with increasing tem- 
perature. Above the coherence temperature, T* ~ 100 K, the 
Drude peak disappears, due to the destruction of Fermi liquid 
quasiparticles. The inset shows the effective number of charge 
carriers iV e ff(u>) per lattice site, calculated from the spectrum 
using Eqns. Q and J5). 



of the optical conductivity ci(w) for several tempera- 
tures. The results are obtained from DMFT using it- 
erative perturbation theory At the lowest temper- 
ature three features are found: a Drude peak at u> = 0, 
a broad absorption band at u> ss 2000 cm -1 ~ U, and 
a band at U/2. The broad band is a result of elec- 
tronic transitions between the Hubbard bands separated 
by U = 10\t 2 \ = 0.3 eV, while the band at U/2 results 
from transitions between the quasiparticle peak and Hub- 
bard bands. 

Single crystals of k-(BEDT-TTF) 2 Cu[N(CN) 2 ]- 
Br^Cli-j; are investigated by polarized reflectivity 
spectroscopy in a broad freq uency and temperature 
range (for details see Ref. Il5l ). The optical response 
can be separated into contributions from correlated 
charge carriers, intradimer charge-transfer transitions, 
and vibronic modes activated by these charge-transfer 
transitions. While the spectra of all the compounds are 
similar at high temperatures, a clear distinction can be 
made below 50 K. A Drude peak is present for a Br 
content of x = 0.73, 0.85, and 0.9 indicating metallic 
behavior. In contrast, no Drude peak in seen in the 
materials with x — and 0.4, consistent with a Mott 
insulating ground state [15]. In order to compare the 
measurements to DMFT calculations, we followed the 
procedure used in Ref. [H to subtract the contributions 



of intradimer electronic transitions (around 3400 cm -1 ) 
and vibrational features from the optical conductivity 
spectra, since the associated physics is not incorporated 
in the model ((T|). 

The optical conductivity, a\{uj), is shown in Fig. [5] for 
the compound with 73% Br. We confine ourselves to 



FIG. 2: (Color online). Frequency dependence of the optical 
conductivity due to the correlated charge carriers of k-(BE- 
DT-TTF) 2 Cu[N(CN) 2 ]Bro.73Clo.27. Experimental data HH 
are shown for conductivity in the layers (E || c) at different 
temperatures. This data compares favourably with theoret- 
ical calculations presented in Fig. [1] The inset shows the 
effective number of charge carriers N c s{u)) per dimer. 



the polarization _£7 1 1 c; [the crystal c-axis is parallel to 
the diagonal k y — k x in the model ©] similar results 
are obtained along the second intralayer crystal axis. At 
T = 150 K, one broad absorption feature is observed 
centered around 2000 cm -1 . When the sample is cooled 
down, the SW shifts towards lower frequencies and be- 
low 50 K a narrow Drude-like component develops. In- 
terestingly, the overall SW decreases for T > 50 K. The 
redistribution of spectral weight with varying tempera- 
ture, observed in Fig. [2] is consistent with the calculations 
shown in Fig. [T] The agreement is particularly impres- 
sive as there are really only two parameters in the calcu- 
lation, the magnitude of i 2 (which sets the vertical scale) 
and the ratio U/t 2 which determines the relative weight 
in the Drude peak. However, we do note two discrep- 
ancies. First, the band at ui ~ U/2 does not appear in 
the experimental data. Second, the temperature scale 
on which the Drude weight collapses is about a factor 
of two larger in the calculation than in the experiment. 
Such discrepancies are not unusual in systems exhibiting 
Kondo physics, as occurs in DMFT MJ. 

To describe the redistribution of spectral weight with 
varying temperature and correlation strength, it is use- 
ful to consider restricted f-sum rule [2, l2l| based on the 
integral: 



MA) 



0i (u>) do; 



2m* 



(3) 



where A is a high-frequency cut-off which excludes in- 
terband transitions (we take A = 3500 cm -1 ), e is the 
electronic charge, n is the total electron density (exclud- 
ing filled bands) and the last identity defines a particular 
effective mass rn* um . Expression ([3]) can be compared to 
the full f-sum rule which is independent of temperature 
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and the strength of the correlations [l|, 0- When the 
conductivity is integrated up to infinity, it can be related 
to the plasma frequency of an electron gas with the total 
electron density. 

In order to make more detailed statements we have to 
consider the effective mass and effective carrier concen- 
tration separately. For a single two-dimensional metallic 
band of non-interacting electrons all the spectral weight 
is in the Drude peak and the sum rule is 



U J OC 



d 2 k d 2 e(k) 
(2irh) 2 dk 2 c 



tine 



. (4) 

where b is the inter layer spacing, and integral [3] is over all 
the occupied states (i.e. the Fermi sea). The last equality 
defines an optical band mass, m^opt [22j |. For the tight- 
binding model @ at half-filling, evaluating Eq. (Tj[]), with 
< 2 = —0.03 eV and t\ — 0.8^2 yields an optical band mass, 
Wfc i0p t = 2.3m e . This agrees very well with rrih.qnt = 
2.5m e obtained from our experimental data (Til. Ila|. 

For a Hubbard model with only nearest-neighbor hop- 
ping on an isotropic lattice this integral is proportional 
to the expectation value of the kinetic-energy operator 
0, H, [5l| . This sum rule even holds in the absence of a 
Drudc response and in the Mott insulating phase. The 
total spectral weight - as calculated in Eq. ((3J - can 
then vary with temperature and the strength of the elec- 
tronic correlations [Bj]. However, we should point out 
that although this connection with the kinetic energy will 
hold qualitatively for our Hubbard model, it will not hold 
quantitatively because the system of interest is spatially 
anisotropic. In the Hubbard model fl} there are next- 
nearest neighbor terms (t\) which are comparable to the 
nearest neighbor terms and frustrate the kinetic energy. 

For the further analysis of the effects of electronic cor- 
relations on the charge-carrier kinetic energy, we rewrite 
Eq. (|3|) and define a frequency-dependent carrier density: 



N eS (ld) 



2TOfc 
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0~\ (id') did' 



(5) 



where at infinite frequencies N c g is equal to the number 
of charge carriers on a dimer (f-sum rule), while devia- 
tion from this value indicate the change in kinetic energy; 
i.e. the influence of electronic correlations. For a weakly 
correlated metal, N e g(id) increases rapidly from zero to 
unity once id becomes larger than the intraband scat- 
tering rate, because almost all the spectral weight is in 
the Drude peak (see the U = 0.06 eV curve in Fig. [3]). 
As U/W is increased, both in theory and experiment 
(Figs. [3] and |4|) , SW is transferred to higher frequencies 
and N e s(id) saturates above frequencies of about id ~ U. 
The absolute value of N c g is substantially suppressed. 
Qualitatively, the reduction in the total spectral weight 
with increasing correlations is due to the diminishing ki- 
netic energy of the charge carriers. 

In the inset of Fig. [TJ N cS (u>) evaluated using DMFT 
with U = 0.3 eV is plotted for different temperatures. 
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FIG. 3: (Color online). Effective charge-carrier number per 
(BEDT-TTF)2 dimer approaching the Mott transition. In the 
left panel the results of a DMFT calculation of N e g(id) (at a 
temperature of 50K) show how the spectral weight is redis- 
tributed as one approaches the transition to the Mott insulat- 
ing phase which occurs at U = 15\tn\ = 0.45 eV. The SW re- 
distribution and saturation of N c ft (u) is in agreement with the 
experimental data for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br !c Cli_ !c 
for x = 0.73, 0.85 (at low temperature) shown in the right 
panel. Typical error bars for the experimental data are shown. 



At large frequencies, the integrated spectrum saturates 
to about 0.6. The calculations are consistent with 
experimental results for k-(BEDT-TTF) 2 Cu[N(CN) 2 ]- 
Bro.73Clo.27, shown in the inset of Fig. [2] and the right 
panel of Fig. [3] With increasing temperature, N c g(uj) 
is suppressed at low frequencies and spectral weight is 
transferred to high id, whereas the total spectral weight 
of the correlated carriers is not conserved because of the 
gradual destruction of quasiparticles with increasing tem- 
perature. 

For a more sophisticated analysis of the mass enhance- 
ment of the correlated carriers, the frequency dependent 
mass renormalization can be quantified with the gener- 
alized Drude model [H, [23| : 



(2Mla(A) 



\/r(id) — iidm* (iv) j 'mb, 



cr(id) = ai(id) + ia 2 (id) = 

j , u; j - /u.7/; ut' i ■ /;/;, ,„>. 

(6) 

where m*(u) and 1/t(u>) are a frequency dependent op- 
tical effective mass and scattering rate, respectively. At 
high frequencies the effective mass becomes equal to the 
optical band mass, since at high frequencies, u(uj) is 
dominated by the imaginary part. The frequency depen- 
dences of this scattering-rate and effective-mass enhance- 
ment are plotted in Fig. SI A sharp drop in m*(uj) and a 
peak in 1/t(co) both occur in DMFT at a frequency, u c 
as shown in Figs. 4(a) and 4(b). Below a frequency scale, 
id* Fermi liquid behavior, i.e., 1/t(io) oc uj 2 , occurs. For 
U = 0.3 eV, the DMFT calculations give: id* « 400 cm" 1 
(inset of Fig. 3^,) in good agreement with the experimen- 
tal value. Both u c and id* are shifted to lower frequencies 
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FIG. 4: (Color online). Frequency dependence of the scat- 
tering rate and effective mass extracted from an extended 
Drude model analysis of the optical conductivity, (a) and (b) 
shows the results of DMFT calculations for different U/W 
and T = 50K. (c) and (d) show the experimental results for 
the same quantities at T = 5K. The inset of (a) shows how 
below a frequency us* ta 400 cm -1 the scattering rate found 
for both experiment and theory has the quadratic frequency 
dependence expected for Fermi liquid quasiparticles. Clearly 
as the Mott insulating phase is approached the effective mass 
and the scattering rate increase significantly. 



with increasing U/W, consistent with the shift observed 
as x decreases from 0.85 to 0.73. The collapse of the effec- 
tive mass enhancement at high frequencies is reminiscent 
of what occurs in heavy fermion compounds Q and the 
'kink' seen in ARPES experiments on cuprate supercon- 
ductors and has recently been argued to be a ubiquitous 



feature of strongly correlated electron systems [24], |25| 



The large increase in the effective optical mass as 
the system comes closer to the Mott insulating phase 
(due to increasing CI content) is consistent with the 
large mass renormalization expected from Brinkman- 
Rice and DMFT pictures This is distinctly different 
from what occurs in doped Mott insulators, such as the 
cuprates, for which the effective mass deduced from an 
extended Drude model analysis is weakly dependent on 
the doping as the Mott insulating phase is approached 
2(|[27j]. For U/t 2 = 2,8, 10, 12, the enhancements in the 
effective mass of the Fermi liquid quasiparticles calcu- 
lated from DMFT are m* /m b , op t = 1.1,2.1,3.1, and 4.1, 
respectively. 

In conclusion, quasiparticles at the verge of localiza- 
tion in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br :c Cli_ :c display a 
sizeable redistribution of spectral weight as the temper- 



ature and x vary. The strong local Coulomb correla- 
tions developing close to the Mott transition at x = 0.7 
lead to spectral-weight transfer from low to high frequen- 
cies. The effective masses extracted from optical data are 
strongly enhanced approaching the Mott transition, in 
contrast to doped Mott insulators. These effects lead to 
the destruction of coherent excitations above a frequency 
scale, uj* . 
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